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ABSTRACT 

We calculate the spectra of inverse Compton (IC) emissions in gamma-ray burst 
(GRB) shocks produced when relativistic ejecta encounters the external interstellar 
medium, assuming a broken power-law approximation to the synchrotron seed 
spectrum. Four IC processes, including the synchrotron self-Compton (SSC) processes 
in GRB forward and reverse shocks, and two combined-IC processes (i.e. scattering of 
reverse shock photons on the electrons in forward shocks and forward shock photons 
on the electrons in reverse shocks), are considered. We find that the SSC emission 
from reverse shocks dominates over other emission processes in energy bands from 
tens of MeV to tens of GeV, for a wide range of shock parameters. This mechanism 
may be responsible for the prompt high energy gamma-rays detected by the Energetic 
Gamma Ray Experiment Telescope (EGRET). At TeV energy bands, however, the 
combined-IC emissions and/or the SSC emission from the forward shocks become 
increasingly dominant for a moderately steep distribution of shocked electrons. 



Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal 



1. Introduction 



The current standard model for gamma-ray bursts (GRBs) and their afterglows is the fireball 
shock model (see Piran 1999 for a recent review). It involves a large amount of isotropic equivalent 
energy release, Eq ~ 10 52-54 ergs, within a few seconds and in a small volume with negligible 
baryonic load, which leads to a fireball that expands ultra-relativistically into the external 
medium. A substantial fraction of the kinetic energy of the baryons is transferred to a non-thermal 
population of relativistic electrons through Fermi acceleration in the shock (e.g. Meszaros & Rees 
1993). The accelerated electrons cool via synchrotron emission and inverse Compton scattering in 
the post-shock magnetic fields and produce the radiation observed in GRBs and their afterglows 
(e.g. Paczyhski & Rhoads 1993; Katz 1994; Sari et al. 1996; Vietri 1997; Waxman 1997a; Wijers 
et al. 1997). The shock could be either internal due to collisions between fireball shells caused 
by outflow variability (Paczyriski & Xu 1994; Rees & Meszaros 1994), or external due to the 
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interaction of the fireball with the surrounding interstellar or wind media ( Meszaros & Rees 1993; 
Dai & Lu 1998; Chevalier & Li 1999). 

When the relativistic ejecta encounters the external medium, a relativistic forward shock 
expands into the external medium and a reverse shock moves into and heats the fireball ejecta. 
The shocked ambient and ejecta materials are in pressure balance and separated by a contact 
discontinuity. The forward shock continuously heats fresh gas and accelerates electrons, producing 
long-term afterglows through the synchrotron emission (e.g. Waxman 1997a,b; Wijers et al. 1997; 
Vietri 1997; Wang et al. 2000a,b; Dai & Lu 1999, 2000; Dai et al. 1999; Huang et al. 1998, 2000; 
Gou et al. 2001). On the other hand, the reverse shock operates only once and after that emission 
from the fireball ejecta are suppressed after the reverse shock crosses the ejecta, since then the 
ejecta expands and cools adiabatically. 

When a reverse shock crosses the shell, the shocked shell and the forward shocked external 
medium carry comparable amount of energy. However, the typical temperature of the shocked 
shell electrons is lower since the particle number density is higher. So, the typical frequency of 
the synchrotron radiation from the shocked shell is considerably lower. A strong prompt optical 
flash (Akerlof et al. 1999) and late time radio flare behavior (Kulkarni et al. 1999), accompanying 
GRB990123, have been attributed to the synchrotron emission process from this reverse shock 
(Sari & Piran 1999; Meszaros & Rees 1999). 

In a recent paper (Wang, Dai & Lu 2001; hereafter WDL), we analytically study the high 
energy 7-ray emission from the SSC process in the reverse shocks. The result showed that this 
emission dominates over the components of the the synchrotron and IC emissions from forward 
shocks at high energy 7-ray bands from tens of Mev to tens of GeV. Furthermore, we suggest that 
this mechanism may explain the observations of the prompt high energy gamma-rays detected by 
EGRET, such as those from GRB930131. 

In this work, we numerically calculate the SSC radiation components in the forward and 
reverse shocks. We also consider another two combined-IC processes, i.e. scatterings of reverse 
shock photons on the forward shocked electrons and forward shock photons on the reversely 
shocked electrons. For a wide range of shock parameters, the present numerical results confirm 
our previous suggestions that the SSC emission from reverse shocks is the most important at the 
energy bands to which EGRET is sensitive. 

While IC emissions from afterglow forward shocks (e.g. Sari et al. 1996; Totani 1998a; 
Waxman 1997; Panaitescu & Meszaros 1998; Wei & Lu 1998, 2000; Chiang k Dermer 1999; 
Dermer, Bottcher &: Chiang 2000; Dermer, Chiang &: Mitman 2000; Panaitescu & Kumar 2000) 
and GRB internal shocks (Papathanassiou & Meszaros 1996; Pilla Sz Loeb 1998; Panaitescu & 
Meszaros 2000) have attracted a lot of investigation, the studies for IC process in reverse shocks 
and combined-IC scatterings between reverse and forward shocks (i.e. very early afterglows) are 
still quite preliminary (see e.g. Meszaros, Rees & Papathanassiou 1994). 

In section 2, we briefly describe the hydrodynamics evolution of a fireball shell. We present 
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in section 3 the synchrotron seed spectra and shocked electron distributions, then perform a 
numerical calculation of IC emissions in section 4. Finally, we give a summary and discussion. 

2. Hydrodynamics of a fireball shell 

Let's consider an ultrarelativistic cold shell with energy E, initial Lorentz factor r] and a width 
A in the observer frame expanding into a cold external interstellar medium (ISM) . When the shell 
sweeps up a large volume of ISM matter, it begins to be decelerated significantly. The interaction 
between the shell and the ISM matter leads to two shocks: a forward shock propagating into the 
ISM and a reverse shock moving back into the shell. There are four regions separated by the 
two shocks and a contact discontinuity: the cold ISM (denoted by the subscript 1), the shocked 
ISM (2), the shocked shell material (3) and the unshocked shell material (4). From the shock 
jump conditions and the equalities of pressure and velocity along the contact discontinuity, we can 
evaluate the Lorentz factor 7, the pressure p and the number density n in the shocked regions in 
terms of three variables ni, and i] (Blandford & McKee 1976; Sari & Piran 1995). 

Sari & Piran (1995) showed that, if the shell is thin: A < Z/2r/ 8 / 3 , where I = (3E , /47rnim p c 2 ) 1 / 3 
is the Sedov length, the reverse shock is Newtonian (with shell spreading, the reverse shock can be 
mildly relativistic) , which means that the Lorentz factor of the shocked shell material 73 is almost 
unity in the frame of the unshocked shell and 73 ~ rj. On the other hand, if the shell is thick: 
A > Z/2r/ 8/3 , the reverse shock is relativistic with 73 ~ 77/273 and it then considerably decelerates 
the shell material (Kobayashi 2000). 

After the shell being decelerated, the evolution of the adiabatic forward shock follows 
the Blandford-McKee solution (Blandford & McKee 1976). For the shocked shell, its dynamic 
evolution may be complex. Kobayashi & Sari (2000) showed that if the fireball shell has relativistic 
temperature, the Blandford-McKee solution can be regarded as an adequate description quite 
early on. But for the case that the reverse shock is only mildly-relativistic, the evolution of the 
shocked shell may deviate the Blandford-McKee solution. This may affect the time evolution 
of the late-time reverse shock emission, but does not affect our computation of the peak-time 
emission. 

3. Synchrotron emission spectra and electron distributions in GRB shocks 

3.1. The synchrotron seed spectrum 

The IC emissions depend on both the seed photon spectrum and the shocked electron 
distribution. We assume the synchrotron seed spectrum to be described by four broken power-law 
segments as given in Sari et al. (1998). The distributions of newly-shocked electrons in both the 
forward and reverse shocks are assumed to be a power law of index p (Nfr) oc 7™ p ), with the 
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minimum Lorentz factor of shocked electrons in the shell rest frame being 7 m = ^ E |rf£e7s/i, 
where £ e is the fraction of thermal energy carried by electrons (Sari, Piran & Narayan 1998). Here 
for the forward shock 7^ = 72 = 73 = V and j s h = 73 for the reverse shock. Assuming that £g 
is the fraction of the thermal energy carried by the magnetic field, the magnetic field strength 
is B' = 12G6j_ 2 1/2 (7 2 /300)n^ 2 , where £ B ,-2 = 6?/10~ 2 , and m, = m/10 . The shocked 
relativistic electrons cool through synchrotron emission and IC scatterings of the synchrotron 
photons; thus the Lorentz factor of the electrons that cool on a timescale equal to the dynamic 
timescale is given by 

_ 67rm e c(l + z) 4500 ^ 72 s_ 3 .-i, 1 + z 1 

^~ (Y + l)a Tl2 {B>)H "y + l^- 2l 300 j ll 2 Jl > ' [ ) 

where t\ = t/lOs denotes the time in the observer frame, Y is the Compton parameter, z is the 
redshift of the burst source and ut is the Thomson scattering cross section. The Compton 
parameter Y, defined to be the ratio between the inverse Compton luminosity of electrons and the 
synchrotron luminosity, can be expressed as Y = Ljc/L syn = U ra d/UB ~ (— 1 + \A + Arj e S, e / ^b) 1^-, 
where U syn and Ub are the energy density of synchrotron radiation and magnetic field , 
respectively, and rj e represents the fraction of the electron energy that was radiated away ( Sari 
& Esin 2001; Panatescu & Kumar 2000). For fast cooling electrons, Y is estimated to \j£, e /£,B 
and it may be smaller for slow cooling electrons. However, as the two combined-IC processes (i.e. 
scattering of reverse shock photons on the electrons in forward shocks and forward shock photons 
on the electrons in reverse shocks) are considered here, the IC cooling of the electrons in both 
shocks should be dominated by the larger of the two photon fields. Because both shocks have 
comparable energy, the magnetic filed energies in the two shocks are comparable, and therefore Y 
in two shocks is almost the same and equals to y/£ e /£,B as long as one shock is in the fast cooling 
regime. 

According to the shell thickness, we divided our discussion into two cases: the thin shell case 
and the thick shell case. The thin shell case implies that the deceleration timescale td ec (defined 
as the observer time at which the heated ISM energy is comparable to the initial energy E), i.e. 

r dec n _l ^ _ +^ F l/3„-l/3„-8/3, 

is larger than the shell crossing time (A/c), where -E53 = £^/10 53 erg and 17300 = ??/300, while the 
thick shell has A/c > td ec - The reverse shock emission peaks at td ec for the thin shell case and at 
A/c for the thick shell. Our following calculations of the IC spectra correspond to this peak time. 



**c = S(l + *) = m^)E^n^%oT b, (2) 



3.1.1. The thin shell case 

As described by Sari et al. (1998), the synchrotron radiation from the shocked electrons can 
be approximated by a broken power-law spectrum with three characteristic break frequencies. One 
is the self-absorption frequency, v a . The other two are the peak frequencies of the emission from 
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the electrons with the characteristic Lorentz factor j m and the cooling Lorentz factor j c , denoted 
as v m and v c , respectively. As usual, the fraction of thermal energy carried by electrons, £ e , and 
magnetic field, (g, is assumed to be similar in the forward and reverse shock. At the peak time 
tdec °f the reverse shock emission, the break frequencies v m from the two shocks are, respectively, 

vZ = 1 T ^ 2eB ' = 6.4 x 10 15 (^) 2 (%^ /2 a^ft— ) Hz, (3) 
m 1 + z 2imi e c K p - V v 0.6 7 SB ~ 2 ,300 1,0 v l + z> ' v; 

and 

jjs _ 1 T (l / L S ) 2 eB' _ 20 p - 2 2 g e 2 1/2 4 1/2, 2 

The equipartition values are chosen to be the ones inferred for GRB970508: £ e = 0.6 and £b = 0.01 
(Granot et al. 1999; Wijers & Galama 1999). Similarly, the cooling break frequencies of both 
shocks are given by 

Vc ~ Vc "1 + z 2nm e c " (Y + 1) 2 53 ^'-^W 4 + z j ^ l&j 
The peak flux of the synchrotron radiation is given by (Wijers & Galama 1999) 

where P' Vm is the peak spectral power and Dl is the luminosity distance of the burst. Note that 
for the forward shock, N e is the total number of swept-up ISM electrons while for the reverse one, 
iV e is the total number of shocked shell electrons. Then, we obtain the peak flux of the two shocks 
at the time tdec- 

C = 4.8 (^)^ 2 28 C 2 Wo-:^5 5 3 /4 ( T ^^)- 3/4 Jy, (7) 
fit = 26 (^)D2%e^ 2 E 5 sn^ mJy, (8) 

where L> L , 28 = D L /10 28 cm. 

Now we derive the synchrotron self-absorption frequency of the reverse shock emission. In the 
comoving frame of the shocked gas (denoted with a prime), the absorption coefficient a' u , scales 
as a' u , oc i/~( p+4 )/ 2 for v 1 > u p = min(f^ ra ,^) and as a' u , oc z/~ 5 / 3 for v' < u p . In this frame, the 
absorption coefficient for v 1 < u p is given by (Rybicki & Lightman 1979) 



K , _ (^f {m ^KB^r (?E±1) r f *±*} f 4) " 5/3 , 

87rm e V27rm3 c 5 y v ; \ 12 J \ 12 J p \v' p j 

(9) 

where K = (p — 1)714(473 + 3)7^~ 1 , 7 P = min(7 m ,7 c ), e is the electron charge and T(x) is the 
Gamma function. Noting that the shocked shell width in the comoving frame Ar' = r/rj{A^j^ + 3) 
and setting 

r^) = a/, Ar' = 1, (10) 
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we get the synchrotron self-absorption frequency v r a s in the observer frame, 

v? = 1.4 x lO^f & 5 (||)- 1 ^ 5 - 2 <o 10 ( T ^)Hz 



(11) 



for p = 2.5. Similarly, one can obtain the synchrotron self- absorption frequency of the forward 
shock emission, 



= 2.2 x l^El' W ^ n ^\^)-V\Y + l)^)^. 



(12) 



Please note that because the forward shock always locates in front of the reverse shock, there will 
be almost no detected reverse shock emission below the self-absorption frequency of the forward 
shock. Thus, the synchrotron spectrum of the slow-cooling reverse shock is described by 
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while for the forward shock, which is in the fast cooling regime for typical shock parameters, 
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3.1.2. The thick shell case 

In a thick shell case, the reverse shock becomes relativistic before it crosses the entire shell 
and begins to decelerate the shell material. It crosses a thick shell at T ~ A/c and the peak time 
of the reverse shock emission is comparable to the GRB duration T. The Lorentz factor of the 
shocked shell scales with time as ( Sari 1997) 

73 = (^) 3/8 (^)" 1/4 , (15) 
thus the Lorentz factor of the reverse shock at this peak time is 

73 = 77/273 = 2r ?5 oo^ 3 1/ ^ 1 ^ 10 Q s ) 3/ yff- ( 16 ) 

The reference values of the comoving width and the initial Lorentz factor of the thick shell have 
been chosen to be A = 3 x 10 12 cm and r\ = 500, respectively. The shock radius at the peak time 
of the reverse shock emission is 

r = 2 73 2 cT/(l + z) = 5 x 10 16 ^^(^^s) 1 / 4 cm. (17) 



- 7- 



Using the expressions of 7 m , j c , Y parameter, and Eq.(7), we obtain the break frequencies and 
peak flux of the reverse shock synchrotron spectrum: 

C = 1-8 x 10 16 (|^|) 2 ( 7 ||) 2 ^ 2 - 27 7|oo-ifo 2 ( T ^) Hz, (18) 



_ 7.6 X 10 3/2 -1/2 -L, T 1/2 2 1/2 



(y + i) 2 s ^>- 2 ^hxv v i + z- 
C = o.ac^wV^^lX 7 ^!^ i^) _3/4 J y- (2°) 

Noting that the shocked shell width in the comoving frame Ar' = A 72/ (473 + 3) and the number 
density 723 = (473 + 3)£ , /(47rm p c 2 r/72Ar 2 ), we can also obtain the self-absorption break frequency 
v a using Eq.(10) for the thick shell case. Since for a thick shell, the reverse shock is generally in 
the fast cooling regime, the synchrotron spectrum is described in a way similar to Eq.(14). 



3.2. Electron distribution in shocked material 

The electron distribution Nfa) in shocked shell or shocked external medium is determined by 
the initial shocked electron distribution ^(7) and by electron cooling effect through synchrotron 
and IC radiation and possibly by the self-absorption of the synchrotron photons. If the 
newly-shocked electrons with typical Lorentz factor 7 m cools faster than the shock dynamics 
timescale, the resulting electron distribution takes the form 




if 7 C < 7 < 7 m _ 
_1 if 7 m < 7 < -y max 



where 7 max is the maximum Lorentz factor of shocked electron, which is determined by equating 
the electron acceleration timescale with the synchrotron cooling timescale (e.g. Meszaros, Laguna 
& Rees 1993), i.e. j max = W 8 B'~ 1 / 2 . In the opposite case, most electrons have a random Lorentz 
factor 7 m , and the electron distribution is 

[I 11 Ic ; ^ Jmax 



4. IC emissions from very early external shocks 

After having obtained the seed photon spectrum and the electron distribution, now we can 
compute the up-scattering emissions of the synchrotron radiation by relativistic electrons in the 
very early external shocks. We consider only the first-order IC and neglect higher order IC 
processes, because a once-scattered typical photon by a typical electron with j e has energy of 
order (hu) com ^ ^ m e c 2 in the rest frame of the second scattering electron. Then we can no longer 
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use Thomson limit to the scattering cross section, and the energy gained of the scattered photon 
in each successive scattering will be reduced due to electron recoil and to the necessity of using 
Klein-Nishina scattering cross section (Sari & Esin 2001). 

For single scattering, the IC volume emissivity in the comoving frame for a distribution N(-y) 
of scattering electrons is given by (Rybicki & Lightman 1979; Sari & Esin 2001) 



where x = v' /Aj 2 ^, f' v , is the incident specific flux at the shock front in the comoving frame, and 
g(x) = 1 + x + 2xln(x) — 2x 2 reflects the angular dependence of the scattering cross section for 
7 e > 1 (Blumenthal & Gould 1970). Noting that f'J c = j'J, c 4Trr 2 Ar'/4:irD 2 and the synchrotron 
flux f' v , = f v , Airr 2 / 'Atk D 2 , where Ar' is the comoving width of the shocked shell or ISM medium 
and D is source distance, we obtain the IC flux in the observer frame 



by transforming Eq.(23) into the observer frame. 

Apart from the SSC scattering processes in the reverse and forward shocks, another two 
combined-IC scattering processes are also present. Because approximately one-half of the photons 
arised in one shock region will diffuse into the another shock region from the point of view of the 
comoving frame, the IC flux Eq.(24) for the combined-IC scatterings should be divided by a factor 
of two. Though the scattered photons move isotropically in the comoving frame, the beaming 
effect makes these photons moving along the direction to the observer. 

Our main calculation results follows: 

i) The IC spectral flux from a relativistic thin shell expanding into an ISM at the deceleration 
time are shown in Fig. 1. Typical shock parameters are used: E = 10 53 erg, £ e = 0.6, £b = 0.01, 
p = 2.5 and n\ = 1. Four IC spectra are displayed in the figure, including two SSC spectra from 
the reverse shock {solid curve) and forward shock (dotted curve), respectively, scatterings of the 
reverse shock photons on the forward shocked electrons (dash-dotted curve) and the forward shock 
photons on the reversely shocked electrons (dashed curve). From Fig.l, it can be clearly seen that 
the SSC from the reverse shock dominates over the other three IC components at gamma-ray 
bands less than a few tens of GeV with a peak around a few MeV. At e ~ 10 — 100 MeV, the SSC 
component of the reverse shock exceeds other IC components by about two orders of magnitude. 
Although the peak location of SSC emission strongly depends on this unknown parameter xi e , 
this emission dominates over other IC emissions for a wide range of shock parameters. 

ii) In Fig.2, we present the energy spectra (vfl c ) of the IC emissions with various shock 
parameters for the thin shell case. We find that a) for a wide range of shock parameters, the 
SSC component from reverse shocks is the most important at energy bands from tens of MeV to 




(23) 




(24) 
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tens of GeV, to which EGRET is sensitive. It dominates over both other IC and the synchrotron 
emission components, b) The SSC spectra from reverse shocks above the turnover (y > v^, ) 
are logarithmically more flattening than that of the seed synchrotron spectrum and varying 
continuously . For small values of p, we even find that the IC energy output peak well above 
min(^ s ' /c ", v™' IC \ which is the peak emission frequency for the approximate power-law IC 
spectrum. A similar result was reached by Sari & Esin (2001) for the SSC of afterglow emission. 
Moreover, for small value of p (e.g. p = 2.2), the SSC emission from the reverse shock dominates 
over the synchrotron and IC processes even in the TeV energy bands (see Fig. 2(d)). c) Fig. 2 
also suggest that strong TeV emission should also be emitted from the two combined-IC and 
forward shock SSC processes for most GRBs. For a moderate steep distribution of the shocked 
electrons (e.g. p = 2.5), the combined-IC and/or forward shock SSC become increasing dominated 
at TeV bands. However, it would only be detected from nearby, low-redshift bursts for which the 
attenuation due to intergalactic infrared emission is small. 

iii) EGRET has detected prompt emission above 30MeV from several bright GRBs triggered 
by BATSE (Catelli et al. 1998), among which GeV photons have been detected from GRB930131 
(Sommer et al. 1994; Ryan et al. 1994) and GRB940217 (Hurley et al. 1994). GRB940217 even 
exhibits delayed GeV emission 90 minutes after the trigger. Several models have been proposed to 
explain the delayed and prompt GeV emissions. For example, Meszaros &: Rees (1994) proposed 
that the impact of a relativistic wind from the central engine on the external matter may cause 
the delayed GeV emission; Katz (1994) suggested that the impact of the fireball on a dense clouds 
could produce high-energy gamma-ray emission via 7r° decay process; Vietri (1997) and Totani 
(1998a, b) suggested that the synchrotron radiation of the protons may be responsible for the GeV 
emissions; etc. In a recent paper (WDL), we have suggested that the SSC emission from the 
reverse shock could explain both the flux level and the spectrum of the high energy gamma-rays 
detected by EGRET. We here compute the slope of the photon spectrum at high energy bands 
and plot it in Fig. 3. We can see that at energy bands from tens of MeV to tens of Gev, the 
photon spectrum index a (the photon number dn(hv)jdv oc v a ) ranges from 1.7 to 2.15, which 
is consistent with the observed high energy gamma-ray photon spectrum by EGRET from some 
bright GRBs (e.g. Sommer et al. 1994). 

iv) Then, we investigate the SSC emission from the reverse shock of a thick shell. We 
compute the time-integrated SSC energy spectrum for a thick shell with typical parameters as 
A = 3 x 10 12 , E = 10 53 erg, p = 2.5 and n\ = 1, but with different values of the initial Lorentz 
factor r/: 7] = 300 (dash- dotted curve), rj = 500(dashed curve) and r] = 1000 (dotted curve). They 
are plotted in Fig. 4. in comparison with the thin shell case (solid curve). As r\ increases, the peak 
frequency of the time-integrated SSC energy spectrum increase accordingly. The peak frequency 
of SSC emission of the thick shell is not far from that of the thin shell case with the same i] 
and typical shock parameters, though the reverse shock of the thick shell is relativistic. This is 
because that this IC emission peaks at v™' (this reverse shock is in the fast cooling regime), 
which is close to v r ^ IC for the thin shell. Fig.4 also shows that the peak flux is much lower than 
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the thin shell case. The reason is that at the peak time of the reverse shock emission, the shell has 
traveled to a larger distance, resulting a lower electron scattering optical depth in the shell. 

5. Summary and Discussion 

The fireball model for GRBs involving an ultrarelativistic fireball ejecta expanding into 
an external medium has two successful predictions: one is the afterglow emission from the 
synchrotron process in the forward shock region (Katz 1994; Meszaros &; Rees 1997), and another 
is the prediction of bright optical emission when a reverse shock is present (Sari & Piran 1999b). 
Subsequent observations of multiwavelength afterglows (e.g. Wijers, Rees & Meszaros 1997) and a 
bright optical flash (Akerlof et al. 1999) gave basic confirmation of this model. In this paper, we 
calculate the SSC and combined-IC emissions from these two shock regions at very early phases. 
For a wide range of shock parameters, the SSC emission from the reverse shock dominates over 
the synchrotron and other IC emissions at energy bands from tens of MeV to tens of GeV, while 
the combined-IC and/or the forward shock SSC emissions become increasingly dominant at TeV 
energy bands. 

We further compute the photon spectrum index a (viz. dn{hv)/dv oc v~ a ) of the reverse 
shock SSC emission, which is a ~ 2.0 for typical shock parameters at energy bands to which 
EGRET is sensitive. Based upon this and the SSC spectral flux level, we suggest that this process 
can provide a plausible explanation for the prompt high energy gamma-rays detected from some 
bright bursts, such as those from GRB930131, GRB910503 and GRB940217 etc. In WDL, we also 
derived the decaying light curves of SSC emission after the reverse shock peak time and found that 
it decays quite rapidly, regardless whether the observed band locates above or below the cooling 
break frequency of the inverse Compton component. The planned Gamma-ray Large Area Space 
Telescope (GLAST) mission will have larger effective area and field of view than EGRET, and so 
will likely be able to monitor the time evolution of the high energy gamma-ray flux. 

At TeV energy bands, the combined-IC emissions and the SSC emission from the forward 
shock become dominant over the SSC from the reverse shock for intermediate values of p (e.g. 
p = 2.5). Nevertheless, for small values of p, the SSC from the reverse shock still dominates 
even in TeV energy bands. Recently the Milagro group reported evidence for TeV emission 
from GRB970417a, one of the 54 BATSE GRBs in the field of view of their detector, Milagrito 
(Atkins et al. 2000). An excess of gamma-rays above background is seen during the durations 
of this burst and the chance probability for detecting such an excess is estimated to be less than 
1.5 x 10~ 3 . Totani (2000) suggested that proton-synchrotron model of GRBs provides a possible 
explanation for these observational results. Our calculation is also consistent with detection of 
TeV emission from GRBs that are near enough to avoid serious attenuation due to intergalactic 
infrared radiation field. 

Two currently popular models for GRBs are the mergers of compact objects (neutron stars 
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or black holes) and the cataclysmic collapse of massive stars. In the former model, compact 
objects are expected to have significant spacial velocities so that their mergers would take place 
at many kiloparsecs outside their birthplaces. Thus, GRBs produced by this model would occur 
in the interstellar medium (ISM) with density n ~ 1cm -3 . However, if a collapsing massive star 
(Woosley 1993; Paczyhski 1998) is the origin of the GRBs, the circumburst medium is the wind 
ejected by the star prior to its collapse, whose density decreases outwards. In the wind circumburst 
medium case, the external medium density is much higher than the ISM case at the deceleration 
length-scale and the magnetic field is much higher accordingly. Therefore, from Eq.(2), we 
know that most of the shocked electrons cool to be sub or trans-relativistic ones (7 C — 1 <^ 1) 
on the deceleration time t^ ec (also see Eq.(14) in Dai & Lu 2001). With such low values of j c , 
the self-absorbed cyclo-synchrotron radiation and multiple Compton scatterings of these cooled 
electrons may construct the low-energy part of the resulting spectrum. At the high-frequency part, 
the spectrum is similar to the fast-cooling spectrum of relativistic electrons, i.e. the spectrum is 
still in the form of f u oc v^ 1 ! 2 extending up to v m and then f v oc v~' p upward, produced by the 
cooling of newly-shocked electrons. The resulting IC spectra in this case is therefore much more 
complicated. Also, we estimate that the IC spectra flux at GeV to TeV bands may be as high as 
that of the ISM case. 
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FIGURE CAPTION 

Figure 1. The spectra of the IC emissions at the reverse shock peak time for typical shock 
parameters: E = 10 53 erg, £ e = 0.6, £b = 0.01, p = 2.5 and an ISM external medium with m = 1. 
The solid and dotted curves represent the SSC emissions from the reverse shock and forward 
shock, respectively. Also plotted are the IC emissions of scatterings of reverse shock photons 
on the forward shock electrons (dash-dotted curve) and forward shock photons on the reversely 
shocked electrons (dashed curve). 

Figure 2. The energy spectra of synchrotron and IC emissions at the reverse shock peak 
time for the ISM circumburst environment case with various shock parameters: &)E = 10 53 erg, 
£ e = 0.6, £b = 0.01, p = 2.5 and n x = 1; h)E = 10 52 erg, £ e = 0.6, £ B = 0.01, p = 2.5 and n x = 1; 
c)E = 10 53 erg, £ e = 0.6, £ B = 10~ 4 , p = 2.5 and m = 1; d)E = 10 53 erg, ^ e = 0.6, £ B = 0.01, 
p = 2.2 and n\ = 1. The thin dash-dotted and dashed curves represent the synchrotron spectra of 
the reverse shock and forward shock, respectively. The four IC spectra are shown by the curves in 
the same way as in Fig. 1. 

Figure 3. The high energy gamma-ray photon spectrum index a of the SSC emission 
from the reverse shock with shock parameters as used in Fig.l. 

Figure 4. Comparison of the time-integrated energy spectrum between the thin shell case 
(solid curve) with shock parameters as used in Fig.l and thick shell case with the same shock 
parameters E = 10 53 erg, A = 3 x 10 12 cm, p = 2.5 and n\ = 1 but with different values of rj : 
r] = 300(dash- dotted curve), rj = 500(dashed curve) and rj = 1000 (dotted curve). 
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